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ABSTRACT 

Shipboard  measurements  of  temperature  and  velocity  fluctuations  were 
performed  to  determine  optical  propagation  properties  of  the  marine 

boundary  layer.  Empirical  expressions  describing  the  temperature  struc- 

2 
ture  parameter,  CT  ,  in  terms  of  the  Richardson  number,  Ri ,  overland  were 

used  to  evaluate  data  obtained  for  open  ocean  conditions.  Profiles  of 

the  dissipation  of  turbulent  kinetic  energy,  e,  with  respect  to  height, 

Z,  and  with  respect  to  a  stability  parameter,  Z/L,  were  examined  for 

open  ocean  conditions. 

In  general,  there  was  little  correlation  between  the  measured  CT 

and  the  stability  parameter.  However,  the  distribution  of  e  with  height 

for  both  the  stable  and  unstable  cases  showed  little  deviation  from  that 

expected  for  near  neutral  conditions. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 11 

II.  THEORETICAL  BACKGROUND - 13 

A.  GENERAL - -  13 

B.  STABILITY  CONSIDERATIONS  20 

III.  DATA  COLLECTION 23 

A.  PLATFORM  AND  LOCATION - 23 

B.  INSTRUMENTATION  AMD  PROCESSING  EQUIPMENT  23 

IV.  ANALYSES  PROCEDURES - 33 

A.  ANALYSES  OF  MEAN  DATA 33 

B.  ANALYSES  OF  FLUCTUATION  DATA --  34 

C.  ANALYSES  OF  VARIANCE  DATA 43 

D.  HOT-WIRE  CALIBRATION 46 

V.  RESULTS - --  49 

VI.  CONCLUSIONS - 59 

LIST  OF  REFERENCES 70 

INITIAL  DISTRIBUTION  LIST  - - 71 


LIST  OF  TABLES 


TABLE  I.    Comparison  of  U*  values  -----  -  —  _______  53 

2 
TABLE  II.   CT  and  Ri  results  from  previous  studies  ------  60 

2 
TABLE  III.  C-j-  and  Ri  results  April  cruise 65 

TABLE  IV.   e  results  April  cruise 69 


LIST  OF  FIGURES 

1.  Effects  of  unstable_and  stable  conditions  on 

height  profiles  ofO--------------------  19 

2.  The  dependence  of  Richardson  number  on  stability  ------  22 

3.  Comparison  of  dimensionless  wind  shear  observa- 
tions with  interpolation  formulae  -------------  22 

4.  Optical  ranges  and  R/V  Acania  position  (normal) 

in  the  Monterey  Bay  --------------------  24 

5.  Sensor  location  on  board  R/V  Acania  ------------  25 

6.  Thornwaite  Associates  cup  anemometer  wind  pro- 
file register  system  --------------------  27 

7.  TSI  model  1210  probes  and  wind  vanes  ------------  29 

8.  Hewlett  Packard  model  HP-2350  temperature 

sensitive  quartz  crystal  probes  --------------  30 

9.  Aspirated  shelter  ---------------------  30 

10.  Dunmore  type  lithium  chloride  humidity  sensor  -------  30 

11.  Mean  wind  profile  ---------------------  35 

12.  Potential  temperature  profiles  ---------------  36 

13.  Virtual  potential  temperature  profiles  -----------  37 

14.  Velocity  spectrum  ---------------------  39 

15.  Temperature  spectrum  --------------------  40 

16.  Spectrum  calibration  plot  -----------------  41 

17.  Sample  strip  chart  section  -----------------  45 

18.  In-situ  calibration  plot  ------------------  48 

19.  Wyngaard  et  al  (1971)  prediction  curve  -----------  50 

20.  Observed  prediction  curve  -----------------  51 


21.  Erroneous  profile     ---------------------55 

22.  Garratt's  results     ---------------------56 

23.  Observed  results  --- 57 

24.  e  versus  logZ - 58 

25.  e  versus  log  Z-----------  —  ----------  58 


LIST  OF  SYMBOLS 

2 

C  Refractive  index  structure  function  parameter 
n 

CT  Temperature  index  structure  function  parameter 

e  Rate  of  dissipation  of  turbulent  kinetic  energy 

f  Temporal  frequency 

g  Acceleration  due  to  gravity 

k  Wave  number 

L  Monin-Obukhov  stability  length 

L  Outer  scale,  lower  limit  of  inertial  subrange 

I  Inner  scale,  upper  limit  of  inertial  subrange 

P  Atmospheric  pressure 

r  Separation  distance 

n  Refraction  index 

Ri  Richardson  number 

S(i)  Spectral  density,  i=f  or  k 

T*  Scaling  temperature 

0  Mean  horizontal  wind  speed 

0*  Friction  velocity 

9  Potential  temperature 

9  Virtual  potential  temperature 

Z  Height 

Z/L  Stability  parameter 


ACKNOWLEDGEMENTS 

I  would  like  to  give  special  recognition  and  thanks  to  Dr.  Kenneth 
L.  Davidson  for  his  expert  guidance,  encouragement  and  support  during 
completion  of  this  study.  Many  thanks  also  go  to  Dr.  Thomas  Houlihan 
for  his  technical  assistance  and  advice.  Much  appreciation  is  extended 
to  both  Dr.  Gordon  Schacher  and  to  Dr.  Chris  Fairall  for  their  invaluable 
help  in  the  analysis  of  the  data  collected  for  this  project.   I  thank 
Mr.  Steve  Rinard  for  his  careful  analysis  of  the  computer  programs 
developed  for  this  thesis. 


10 


I.   INTRODUCTION 

The  advent  of  more  complex  and  complete  numerical  models  for  analyz- 
ing and  predicting  large  scale  atmospheric  motion  has  also  led  to  better 
specifications  of  the  boundary  layer.  The  boundary  layer  is  that  region 
in  the  first  kilometer  over  the  sea  which  is  defined  principally  by 
turbulent  transfer  of  momentum  and  heat. 

Optical  propagation  through  the  atmosphere  is  affected  by  the  re- 
fractive nature  of  the  medium.  In  addition  to  the  regular  variation  of 
atmospheric  refractive  index  with  height,  there  exist  small  inhomogenei- 
ties  in  the  refractive  index  associated  with  fluctuations  in  the  tempera- 
ture and  velocity  of  the  air.  These  fluctuations  cause  random  phase 
and  amplitude  distortions  in  propagating  wave  fronts  and  thus  degrade 
spatial  and  temporal  coherence  in  the  transmission.  The  magnitude  of 
the  effects  places  limitations  on  optical  system  performance  and  must  be 
included  in  design  and  execution  considerations. 

Descriptions  of  the  small  scale  fluctuations  which  affect  optical 
propagation  have  not  been  as  complete  nor  in  the  quantity  for  the  over- 
water  regime  as  for  the  overland  regime. 

Initial  experimental  efforts  to  verify  turbulence  theory  predictions 
were  conducted  over  land.  Measurements  of  wind  speed,  temperature,  and 
humidity  in  those  investigations  were  made  from  a  stable  platform  with 
auxiliary  instrumentation  well  protected  from  external  influences. 
Under  these  conditions  relatively  accurate  measurements  of  wind  speed, 
temperature  and  humidity  were  performed  under  a  variety  of  weather  condi- 
tions. 


11 


The  marine  environment  on  the  other  hand  presents  a  number  of  prob- 
lems, especially  in  accessibility,  platform  stability,  and  sensor  pro- 
tection. Overwater  descriptions  are  necessary,  even  though  considerable 
progress  has  been  made  in  overland  investigations.  The  necessity  exists 
because  of  the  increasing  evidence  of  the  influence  on  atmospheric 
motions  by  features  unique  to  the  oceanic  environment. 

It  is  the  purpose  of  this  paper  to  compare  overland  and  overwater 
results,  also,  to  evaluate  the  assumption  that  near  neutral  conditions 
apply  for  both  stable  and  unstable  conditions  when  evaluating  the  varia- 
tion of  turbulent  kinetic  energy  dissipation  with  height. 
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II.  THEORETICAL  BACKGROUND 

A.  GENERAL 

Turbulence  properties  of  interest  in  optical  propagation  are  those 
which  could  be  used  to  describe  the  intensity  of  fluctuations  in  the 
refractive  index  and  that  which  could  be  used  to  describe  the  size  or 
scale  of  the  refractive  index  inhomogeneities.  The  former  has  been  re- 
lated empirically  to  scintillation,  beam  wander  and  beam  spread,  and 
the  latter  along  with  the  former  to  the  image  resolution  which  is  im- 
portant with  respect  to  seeing  conditions.  Similarity  expressions  re- 
lating mean  conditions  to  properties  of  the  turbulent  regimes  and  hence 
the  temperature  structure  function  and  the  dissipation  rate  of  turbulent 
kinetic  energy  have  evolved  through  a  series  of  hypotheses  and  asso- 
ciated measurements. 

An  important  consideration  within  turbulence  theory  is  whether  the 
turbulence  is  anisotropic  or  isotropic.  In  general,  anisotropy  exists 
within  large  scale  eddies  and  isotropy  exists  within  small  scale  eddies. 
The  universal  formulae  considered  will  be  valid  only  for  isotropic 
turbulence. 

Anisotropic  turbulence  is  yery   empirical.  Every  time  the  boundary 
changes,  the  turbulence  changes.  However  at  high  wave  numbers,  which 
corresponds  to  small  scale  eddies,  the  turbulence  should  be  independent 
of  the  boundaries.  Kolmogorov  (1941)  postulated  that  at  high  wave 
numbers  turbulence  should  be  isotropic  even  though  it  is  generated 
from  and  embedded  in  anisotropic  eddies  of  low  wave  numbers.  If 
Reynolds  numbers  are  high  enough,  turbulence  will  adjust  through  inertia! 
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transfer  and  viscous  dissipation  until   a  statistical   equilibrium,   inde- 
pendent of  initial   conditions,   is  attained.     This  equilibrium  consists 
of  energy  transfer  into  the  high  wave  number  range  that  is  equal   to 
energy  lost  through  viscous  dissipation. 

On  the  basis  of  the  isotropic  nature  of  small   scale  fluctuations, 
only  one  parameter  is  necessary  to  describe  the  intensity  of  the  refrac- 
tive index  fluctuation  over  many  scales.     It  is  the  refractive  index 

2 
structure  function  parameter,  C     ,  defined  as 


C  2  =   [<n(x)  -  n(x+r)>]2/r2/3  (1) 


where     n(x)     and     n(x+r)     are  refractive  indices  at  two  points  on  a 
line  oriented  normal   to  the  mean  wind  direction  and  separated  by  the 
distance     r    which  is   less  than  the  outer  scale,   L   ,   the  lower  end 
of  the  inertial   subrange,  and  greater  than  the  inner  scale,  I  ,   the 
smallest  scale  of  naturally  occurring  turbulence. 

A  parallel   expression,  which  defines  the  temperature  structure 
function  parameter  is 


CT2  =   [<T'(x)   -  T'(x+r)>]2/r2/3  (2) 


where     T'(x)     and     T'(x+r)     are  temperature  fluctuations  at  two 
points  separated  by  the  distance     r  .     The  refractive  index  is  deter- 
minted  primarily  by  density  fluctuations  and  can,   therefore,  be 

related  to  temperature  fluctuations,   neglecting  himidity. 

2  2 

C         is  related  to     CT     by  the  following 


Cn2  =   [79.xl0"6  P/T2]2  CT2  (3) 
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2        2 
Fortunately,  both  C    and  CT   are  readily  measurable  by  optical 

and  meteorological  means,  respectively.' n  alternate  relationship  for 

2 
CT  ,  which  involves  measurement  of  the  rates  of  dissipation  of  turbu- 
lent kinetic  energy,  e,  and  temperature  variance,  x>  1S 

CT2  =  6Xe"1/3  (4) 

where  3  is  an  empirical  constant  with  a  value  of  3.2.  This  last  form 

2 
enables  indirect  estimates  of  CT   to  be  made  from  mean  conditions 

since  e  and  x  are  easily  related  to  boundary  layer  fluxes  and, 

eventually,  profiles  if  steady  horizontally  homogeneous  conditions  exist, 

2 
Expressions  which  relate  CT   to  mean  properties  of  the  boundary 

layer  are  desirable  because  the  small  scale  measurements  are  impractical 
to  obtain  in  most  operational  or  tactical  regimes  and  situations. 

Since  turbulence  is  nearly  synonymous  with  temperature  fluctuations, 
it  is  ultimately  desirable  to  describe  mean  thermal  stratification  in 
terms  of  atmospheric  bulk  stability  parameters  such  as  the  Monin-Obukhov 
length  (L).  In  this  regard,  measurements  of  both  atmospheric  mean 
profiles  and  fluxes  are  required  for  a  complete  determination  of  atmos- 
pheric transmission  behavior.  Profile  (86  /8z  and  3U/az)  and  bound- 
ary flux  (U*  and  T*)  parameters  appear  in  the  following  expressions  for 
Ri  and  L 


gae/az 

Ri  =  1 _  (5a) 

9(3U/8z) 


L  =  SU*2/gKT*„  (5b) 
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The  following  similarity  predictions  for  the  dependence  of  e  and 

X  on  momentum  and  heat  fluxes  and  height  were  considered  by  Wyngaard, 

2 
et  al.  (1971)  in  deriving  an  empirical  expression  for  estimating  CT 

from  mean  stability  parameters 


cZ/U*3  =  f^Z/L)  (6a) 


XZ/T*U*=  f2(Z/L)  (6b) 


2 
The  form  of  the  empirical  expression  for  Cj   is  obtained  by  direct 

substitution  of  equations  (6a  and  b)  into  equation  (4)  yielding 


CT2  =  T*2Z"2/3f3(Z/L)  (7) 


Furthermore  since  (Z/L)   and  Ri   can  be  functionally  related, 
Businger,  et  al.    (1971),  a  parallel   dependence  on  Ri   can  be  obtained, 
viz. , 


CT2  =  Z4/30e/3z)2<j>3(Ri)  (8) 


The  functions     f,    ,  f „  ,  f ,  ,  and     <J>3     in  equations   (6),    (7),  and 
(8)   are  empirical   and  formulated  on  the  basis  of  observations  of  both 
temperature  fluctuations,  momentum  and  heat  fluxes,  and  mean  gradients 
of  wind  speed,  temperature,   and  humidity. 

The  relation  expressed  by  equation   (8)   provides  a  desired  depen- 

2 
dence  of     CT      on  more  readily  measured  mean  stability  (Z,   89/ 3z,  and 

Ri ) .     The  forms  of     f~     and     4>3     and  the  data  obtained  from  the  exten- 
sive AFCRL  study  of  turbulence  structure  over  a  flat,  unobstructed 
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Kansas  plain  are  presented  in  Wyngaard.  As  will  be  shown  later,  avail- 
able marine  data  does  not  appear  to  agree  for  f3  as  well  as  expected 
with  the  overland  predictions,  i.e.  Hughes  (1976).  The  overland  pre- 
dictions for  f,  and  f~  in  equation  (6),  were  as  follows 

f3(Z/L)  =  4.9(1-7(Z/L))"2/3,  Z/L  <  0 

f3(Z/L)  =  4.9(1+2.75(Z/L))  ,  Z/L  >  0       (9) 

Small  scale  velocity  fluctuation  properties  are  of  interest  in 
optical  propagation  because  image  resolution  has  been  empirically  re- 
lated to  the  innerscale,  £  ,  which  is  defined 


lo   =  ^/£)V4  (]°) 


where    y     is  the  kinematic  molecular  viscosity  and     e     is  the  dissipa- 

2 
tion  rate  of  turbulent  kinetic  energy  similar  to     C,     .     £     can  be  ob- 
tained from  either  one-dimensional  velocity  variance  spectral   estimates 
in  the  inertial   subrange  or  from  velocity  structure  function  estimates. 

e     and     I       can  be  functionally  related  to  mean  profile  and  flux 
estimates   (U*,  T*  and  Ri   or  Z/L)   on  the  basis  of  the  empirical   expres- 
sion presented  previously  (equation  6a).     For  the  purpose  of  examining 
overwater     e     results  from  different  stability  conditions  described  by 
Richardson  numbers   (Ri),  equation   (6a)   can  be  expressed  as  follows 


£  =  U*3/Kz  ^(Ri)  (11) 


Since  Z/L  and  Ri  are  functionally  related. 
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Overwater  e  values  can  also  be  evaluated  on  the  basis  of  U* 
values  estimated  from  them.  Such  estimates  can  be  evaluated  by  compar- 
ing them  with  U*  estimates  computed  from  mean  wind  profile  measure- 
ments. Both  U*  estimates  utilize  the  following  expression  for  the 
mean  wind  gradient 

dU/dz   =  uykz  ^(Ri)  (12) 

where  <j>-,(Ri)  is  defined  in  equation  (11).  <J>-,(Ri)  is  equal  to  1.0 
under  neutral  conditions.  U*  is  assumed  to  be  a  constant  with  height 
in  this  expression  so  it  is  restricted  to  the  constant  surface  layer. 
Integration  of  equation  (12)  for  near  neutral  conditions  yields 


D  =  u*/k  in   Z/ZQ  (13) 


where  U   is  assumed  equal  to  zero,  Z   is  the  roughness  parameter. 
The  logarithmic  profile  suggested  by  equation  (13)  is  altered  by 
stability,  shown  in  Figure  1.  Z   can  be  eliminated  from  equation  (13) 
by  selecting  mean  winds  at  two  appropriate  levels  so  that 


u*  =  k(U2  -  U1)/£n(Z2/Z1)  (14) 

It  is  important  to  note  that  this  expression  relating  U*  to  mean  wind 
values  (Op  and  D-, )  at  two  levels  (Z2  and  Z-. )  is  applicable  only  for 
near  neutral  conditions  (Ri  =  0). 

In  near  neutral  conditions,  turbulent  kinetic  energy  production  is 
assumed  to  be  equal  to  the  rate  of  molecular  dissipation  of  turbulent 
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kinetic  energy  and  the  following  relation  is 
valid 


e  =  U*2  (3D/3z)  (15) 


combining  equations   (12)   and   (15),  assuming  neutral   conditions     4>-.(Ri) 
=  1.0,  and  solving  for    U*  ,  yields 

U*  =   UkZ)1/3  (16) 

Thus,  under  neutral   conditions,  the  friction  velocity  (U*)  can  be  esti- 
mated from  either  mean  wind  profiles  using  equation   (14)   or  from 
fluctuation  data   (involving  turbulent  energy  dissipation)   using  equa- 
tion  (16). 

From  determined  values  of     U*  ,   it  is  also  possible  to  calculate  a 
momentum  drag  coefficient,  C   ,  corresponding  to  a  given  height  in  the 
surface  layer 


Cz  =  U*2/02  (17) 


Numerous  studies  have  been  conducted  to  determine  a  representative 
value  of  C   for  10  meters,  Cardone  (1963).  Check  calculations  of 
C,q  employing  U*  determined  from  both  mean  wind  profiles  and  dissi- 
pation rates  are  performed  in  this  study. 

B.  STABILITY  CONSIDERATIONS 

Observational  experiments  by  Businger  et  al .  (1971)  yielded  a 
definite  relationship  between  the  Richardson  number,  Ri ,  equation  (5a), 
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and  the  Monin-Obukhov  length,  L,  equation  (5b).  Figure  2  from  Businger, 
et  al.  (1971)  lillustrates  this  relationship  quite  well.  The  following 
expressions  are  approximations  of  the  relationship  between  Z/L  and  Ri , 
proposed  by  Dyer  and  Hicks  (1970)  and  Webb  (1970)  for  unstable  and 
stable  conditions,  respectively, 

Z/L  =  Ri  (18) 

Z/L  =  Ri/l-aRi  (19) 

where  a  is  an  empirically  derived  constant  equal  to  5.0. 

It  is  important  to  note  that  Businger,  et  al .  (1971)  observed  the 
Richardson  number  to  approach  a  critical  value  of  0.21  as  Z/L 
approached  °°  suggesting  that  as  stability  increases,  the  flow 
becomes  essentially  non-turbulent.  This  implies  that  for  the  surface 
layer  the  effect  of  mechanical  turbulence  becomes  negligible. 

The  effect  of  hydrostatic  stability  on  surface  layer  mechanical 
turbulence  is  further  illustrated  in  Figure  3  from  Businger,  et  al. 
(1971).  Recalling  equation  (12)  in  light  of  Figure  3,  it  becomes 
apparent  that  the  function  4>m(Z/L)  vaires  relatively  little  with 
respect  to  Z/L  for  unstable  conditions  approaching  a  value  of  0.5. 
However,  the  rapid  increase  of  <f> m(Z/L)  as  the  atmosphere  becomes 
hydrostatically  stable  indicates  that  the  near  neutral  assumptions 
are  invoked  under  only  slightly  stable  conditions. 
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Figure  2.  The  dependence 
of  Richardson  number  on 
stability. 
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Figure  3.  Comparison  of 
dimensionless  wind  shear 
observations  with  inter- 
polation formulae. 
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III.      DATA  COLLECTION 

A.  PLATFORM  AND  LOCATION 

All  observations  were  made  aboard  the  R/V  Acania  while  anchored  in 
about  30  fathoms  of  water  off  Pt.  Pinos  in  Monterey  Bay.  Open  ocean 
differs  from  land  in  the  effects  of  wave  action  on  turbulence,  in  the 
nature  of  aerosols  and  fog,  and  humidity  fluctuations.  These  conditions 
could  best  be  obtained  far  at  sea;  however  the  cost  of  such  activities 
made  it  desirable  to  work  near  land  and  the  optical  propagation  experi- 
ments could  be  coincident.  Pt.  Pinos  and  other  locations  along  the 
Monterey  Bay  shoreline  provided  yery   nearly  the  ideal  situations  for 
both  the  optical  and  meteorological  experiments,  Figure  (4). 

B.  INSTRUMENTATION  AND  PROCESSING  EQUIPMENT 

Measurements  were  made  at  four  levels.  Figure  5  shows  positions 
of  the  towers  which  support  the  instrumentation,  length  of  cable  runs, 
and  position  of  processing  and  recording  equipment. 

The  mean  wind  measurements  were  obtained  with  Thornwaite  Associates 
cup  anemometer  wind  profile  register  systems,  model  number  104.  Three 
cups,  plastic  cones  reinforced  with  aluminum  frames,  are  mounted  on 
stainless  steel  tubes  attached  to  the  main  shaft  at  120°  intervals.  When 
the  cups  rotate  due  to  the  wind  a  slotted  shaft  in  the  anemometer  serves 
as  a  shutter  between  a  light  source  and  a  photoelectric  cell.  The  three 
cup  anemometers  have  the  characteristics  of  low  starting  speeds  with  a 
small  amount  of  internal  friction  which  aids  in  checking  inertial  over- 
shoot. 
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Figure  4.  Optical  ranges  and  R/V  Acania  position  (normal)  in  Monterey 
Bay. 
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Velocity  fluctuation  measurements  were  obtained  using  TSI  model  1210 
probes  with  sliding  support  shields  and  tungsten  wires.  The  shields  per- 
mitted isolation  of  the  sensing  area  for  the  determination  of  the  undis- 
turbed velocity,  V  ,  voltage  reading  before  and  after  each  experiment. 
The  platinum  coated  tungsten  wire  was  small  enough  to  resolve  the  viscous 
dissipation  scale  without  making  corrections  for  wire  length.  The  wire 
had  plated  ends  for  isolating  the  sensing  area  and  thereby  minimizing 
flow  disturbance.  Figure  (6) 

Electronics  associated  with  each  probe  were  a  TSI  model  1054B 
linearized  anemometer  and  a  TSI  model  1056  variable  decade  module.  The 
anemometer  had  a  linear  frequency  response  from  (dead  calm)  0  KHz  to 
10  KHz  and  the  variable  decade  module  operated  with  a  0-60  ohm  range. 

Temperature  fluctuations  were  measured  using  similar  size  sensors. 
A  platinum  wire  was  used  in  place  of  a  tungsten  wire.  The  temperature 
fluctuation  system  was  designed  for  a  resolution  of  0.001 °C  at  frequen- 
cies up  to  1  KHz. 

The  temperature  fluctuations  were  measured  using  a  bridge  developed 
by  GTE  Sylvania,  the  GTE  Sylvania  model  140.  The  system  was  slightly 
modified  for  use  aboard  the  R/V  Acania. 

The  baseband  portion  of  this  system  is  basically  a  balanced  wheat- 
stone  bridge  excited  by  a  3  KHz  signal  with  a  synchronous  detector  on 
the  output.  Segments  of  a  very  small  diameter  platinum  wire  serve  as 
temperature  sensors  in  opposite  arms  of  the  bridge.  The  resistance 
temperature  coefficients  result  in  an  output  from  the  bridge  which  is 

proportional  to  the  temperature  difference  between  the  two  probes.  The 

-4 
sensor  wire  is  0.5  centimeters  long  and  2.5x10   centimeters  in  diameter, 

This  extremely  small  mass  allows  a  response  to  temperature  variations 
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Figure  6.  Thornwaite  Associates  cup  anemometer  wind 
profile  register  system. 
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of  up  to  1  KHz,  while  electronic  amplification  allows  temperature  differ- 
ences as  small  as  0.004°C  to  be  observed. 

Both  wind  and  temperature  fluctuation  data  were  recorded  on  a  San- 
born model  3950  fourteen  channel  tape  recorder.  Real  time  readout  on 
an  eight  channel  chart  recorder,  brush  model  240,  was  used  to  check  the 
quality  of  the  signals  coming  from  the  sensors.  The  chart  data  was  also 
utilized  in  temperature  variance  analyses.. 

Hewlett  Packard  model  HP-2350  temperature  sensitive  quartz  crystal 
probes  were  used  to  obtain  mean  temperature  at  the  sea  surface  and  four 
tower  levels.  Figure  (7)  The  RF  signal  from  the  probes  and  a  reference 
oscillator  were  mixed  in  a  Hewlett  Packard  HP-2801A  readout  unit  to  pro- 
duce a  beat  frequency  whose  signature  can  be  analyzed  to  within  0.001 
°C/Hz.  Each  sensor  received  pre-experiment  calibration  against  a 
platinum  resistance  wire  thermometer  in  a  temperature  controlled  circulat- 
ing water  bath.  It  was  checked  over  the  range  of  expected  temperatures. 
The  accuracy  in  achieving  a  0.005°C  correction  factor  was  a  constant 
for  each  probe.  The  tower  mounted  sensors  were  housed  in  aspirated 
shelters  to  eliminate  radiation  effects.  Figure  (8) 

Mean  relative  humidity  information  was  obtained  using  Dunmore  type 
lithium  chloride  sensors.  Figure  (9)  This  sensor  was  also  placed  in 
the  aspirated  shelter.  The  basic  principle  of  operation  of  this  sensor 
is  resistance  change  in  an  electrolytic  solution  generating  a  reference 
voltage  variation  proportional  to  relative  humidity  change.  Automatic 
temperature  compensation  in  the  instruments  meet  the  following  specifi- 
cations for  relative  humidity,  +  3%  humidity  below  90%  and  +  4%  humidity 
above  90%.  Sensor  calibration  was  accomplished  by  a  comparative  method 
using  a  saturated  saline  solution  in  a  closed  container. 
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Figure  7.  TSI  model  1210  probes  and  wind  vanes, 
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Figure  9.  Aspirated  shelter. 


Figure  10.  Dunmore  type 
lithium  chloride  humidity 
sensor. 
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Figure  8.  Hewlett  Packard  model  HP-2850  temperature  sensitive 
quartz  crystal  probes. 
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Sensor  placement  required  exceptionally  long  cable  runs.  Here  de- 
crease in  system  frequency  response  due  to  their  length  had  little 
effect  in  the  frequency  band  of  interest.  Adjustments  were  made  in 
the  bridges  for  resistance  and  capacitance  of  the  wire  length. 

Data  logging  during  the  experiments  was  accomplished  using  the  NPS 
developed  MIDAS  (Microprogrammable  Integrated  Data  Acquisition  System). 
This  micro-processor  based  data  acquisition  system  utilizes  an  Intel 
8008  central  processor  to  control  the  sampling,  averaging,  and  record- 
ing of  mean  meteorological  data.  All  software  is  written  in  PL/M  to 
facilitate  the  writing  of  a  self-documenting  program. 

The  operator  is  interfaced  with  the  system  via  teletype  for  full 
duplex  input/output  communications  and  program  control  over  the  sample 
start  time  and  the  number  of  samples  to  be  averaged  before  outputting. 
The  operator  may  also  alter  the  present  preset  sample  list  by  adding 
or  deleting  various  sensors  as  they  come  on  line  or  become  inoperable 
or  inconsistent.  Once  initiated,  the  system  is  fully  automated  to 
sample  the  tailored  list  of  sensors  ewery   30  seconds  and  periodically 
print  output  values  averaged  over  the  selected  interval  of  from  one 
minute  to  one  hour.  In  this  study  a  10  minute  averaging  interval  was 
selected. 

Output  values  were  printed  out  on  the  teletype  in  a  columnized  for- 
mat with  the  time  of  print  as  a  leader.  The  teletype  has  a  paper  punch 
incorporated  which  may  be  activated  by  the  operator  to  produce  a  data 
copy  concurrent  with  the  printout.  This  paper  punch  was  also  utilized 
to  produce  data  cards  for  use  with  the  IBM  360  digital  computer.  A 
magnetic  cassette  tape  recorder  has  been  integrated  into  the  system 
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as  a  third  data  output  device.  This  cassette  can  be  interfaced  to  an 

HP  9830  portable  computer  so  that  profile  and  gradient  flux  estimates 

can  be  performed  on  board  automatically  using  the  basic  programming 
capability  of  the  HP  9830. 
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IV.  ANALYSES  PROCEDURES 

A.  ANALYSES  OF  MEAN  DATA 

Mean  wind  data  for  ten  minute  intervals  were  available  from  the 
Naval  Postgraduate  School  developed  MIDAS  (Microprogrammable  Integrated 
Data  Acquisition  System)  output.  MIDAS  read  values  at  ten  second  in- 
tervals and  averaged  them  for  a  ten  minute  period.  Analyses  were  per- 
formed on  data  recorded  on  teletype  printout  and  also  on  the  cassette 
tape. 

The  data  was  edited  for  gross  errors  or  inconsistencies  due  to 
known  instrumentation  malfunctions.  The  criteria  for  retaining  or 
discarding  data  from  individual  levels  or  for  entire  ten  minute  inter- 
val depended  on  performance  check  results  obtained  during  measurements, 
obvious  inconsistencies  between  levels  and  sequential  times,  and  the 
availability  of  fluctuation  data.  After  mean  data  were  obtained  for 
ten  minute  intervals  and  passed  preliminary  editing  for  obvious  erro- 
neous values,  it  was  processed  for  three  applications,  i.e.,  mean  wind 
D,  mean  temperature  T,  and  mean  humidity  q.  These  values  were  plotted 
on  2-cycle  logarithmic  paper.  Since  0,  T,  and  q  are  expected  to  vary 
logarithmically  with  height,  a  best  fit  straight  line  should  describe 
the  distribution  of  data  points.   Best  fit  lines  were  used  to  estimate 
gradient  parameters  and  the  Richardson  number  (equation  5a,  Figures  10 
and  12).  In  general,  the  procedures  were  subjective  and  in  some  in- 
stances, best  fit  lines  with  different  slopes  could  be  visualized  for 
a  given  set  of  data  points.  Therefore,  one  criterion  used  was  not  to 
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give  a  single  data  point  too  much  influence  in  determining  the  best  fit 
line.  Consequently,  the  line  drawn  represented  a  most  probable  posi- 
tion between  data  points  without  necessarily  coinciding  with  any  point 
(Figure  12). 

B.  ANALYSES  OF  FLUCTUATION  DATA 

2 
Spectral  analyses  to  obtain  CT   and  e  values  were  performed  on 

fluctuating  velocity  and  temperature  data  obtained  with  single  wires 
oriented  in  the  vertical  normal  to  the  mean  wind.  These  data  represent 
a  point  measurement  and  yield  temporal  or  time  descriptions  of  the 
fluctuations.  Twenty-one  minute  segments  of  data  originally  recorded 
on  magnetic  tape  were  recorded  into  the  EMR-Schl umber  model  1510  digital 
spectrum  analyzer.  Procedures  for  converting  spectral  values,  obtained 
with  the  analog  spectral  analyzer,  to  engineering  units  and  for  obtain- 
ing turbulence  parameters  from  the  spectra  are  described  in  the  follow- 
ing paragraphs. 

A  necessary  procedure  was  to  scale  the  spectral  plots  in  order  to 
relate  rms  input  voltages  to  power  spectral  densities  (PSD),  variance 
per  unit  frequency.  To  obtain  PSD  levels,  corresponding  to  rms  voltage 
input,  calibrated  scales  charts  had  to  be  constructed. 

For  purposes  of  the  X-Y  plot  format  of  the  analyzer  output,  the  rms 
voltages  were  converted  to  y  =  log-.Q  (voltage)  units  and  a  graduated 
scale  was  constructed  enabling  the  logarithm  of  rms  volts  to  be  inter- 
polated from  the  spectral  plots.  The  vertical  scale,  Y,  values  from 
the  plots  were  adjusted  for  each  spectrum  as  a  function  of  input  and 
spectral  gains.  The  values  were  then  converted  to  PSD  levels  for  use 
in  calculating  turbulence  parameters.  The  relation  used  was 

(cal.  level  V   )2     Y  « 
S(f)  (PSD  units}  =   1t5bandw1dtT   x  (10V   (20) 
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Figure  11.     Mean  wind  profile. 
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Figure  12.     Potential   temperature  profiles. 
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Figure  13.  Virtual  potential  temperature  profiles 
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where  bandwidth  is  the  frequency  range  divided  by  the  number  of 
channels 


Bu  =  Freq  range 

Number  channels 


=  256  Hz/256  =  1  Hz 

The  calibration  level  volts  rms  equals  the  voltage  at  Y=0.  (1  volt  rms 
for  a  3.16  volt  input  setting). 

Amplitude  scaling  calibrations  were  accomplished  using  externally 
generated  white  noise.  Signals  with  a  1  volt  rms  over  a  frequency 
range  of  0.1  Hz  to  1000  Hz,  giving  a  PSD  of  10"3  v2/Hz.  Setting  a 
3.16v  (10  dB)  input  gain  on  the  spectrum  analyzer  insures  that  1  volt 
input  corresponds  to  Y=0.  An  example  of  such  a  calibration  plot  is 

shown  in  Figure  (15). 

2 
The  turbulence  parameters  C.   and  e  were  obtained  from  the 

variance  spectra  on  the  basis  of  the  universal  formulae  for  the  inertia! 

subrange  in  wave  number  space.  These  expressions  predict  a  -5/3  slope 

for  variance  spectra  of  both  velocity  and  temperature  when  plotted  in 

log-log  format.  Figures  14  and  15  are  typical  spectra  considered  in 

analyses.  The  velocity  variance  spectra  have  consistent  -5/3  slope. 

However,  temperature  spectra  often  exhibited  slopes  slightly  different 

from  the  expected  value  of  -5/3.  This  feature  of  temperature  spectra 

has  been  observed  by  others  and  the  existence  of  "cold  spikes"  in  the 

temperature  traces  has  been  postulated  as  a  probable  cause,  Friehe 

(1976). 

Assuming  -5/3  slopes  for  the  variance  specta,  the  intercept  of  the 

slope  with  the  1  Hz  frequency  line  was  the  power  spectral  density  value 

2 
used  in  computing  the  parameter  of  interest,  e     and  C.   . 
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Figure  14.  Velocity  spectrum. 
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Figure  15.  Temperature  spectrum. 
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Figure  16.  Spectrum  calibration  plot. 
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The  measured  PSD  value  obtained  from  velocity  spectra  was  converted 
to  a  spectral  density  in  engineering  units  by  the  equation 


2  6 

Su(f)  =  V  •  PSD 


=  (cm/s/volt)2  •  volt2/Hz  =  (cm/s)2/Hz 


(21) 


where  C  is  the  hot-wire  calibration  factor.  Similarly,  the  equation 
used  to  convert  PSD  values  obtained  from  temperature  spectra  to 
dimensional  values  was 


S  (f)  =  CH2  •  PSD 

=  (°C/volt)2  •  volt2/Hz  =  °C2/Hz      (22) 


where  CH  is  the  calibration  factor  for  the  temperature  system. 

Since  velocity  and  temperature  fluctuations  were  measured  at  a  fixed 
point  in  the  flow,  the  resultant  spectral  estimates  correspond  to 
"temporal"  frequencies,  f,  as  indicated  by  the  notation  S  (f)  and 
Sj(f)  in  the  two  preceding  equations.  To  obtain  e  and  C.   values, 
temporal  (f)  and  space  (k)  scales  must  be  related  in  order  to  use 
equations  (21)  and  (22).  This  is  accomplished  by  Taylor's  (1938) 
"frozen  turbulence"  hypothesis,  that  is 

k  =  2Trf/D  (23) 

where  D  is  the  mean  wind  at  the  measurement  level.  The  term  "frozen" 
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implies  the  turbulence  pattern  remains  unchanged  as  it  sweeps  past 
the  sensor  probe. 

The  following  identies  are  those  used  to  relate  temporal   spectral 
estimates  to  spatial   spectral  estimates 


f-Su(f)   e  k.Su(k)   =  C1   £2/3  k"2/3  (24) 

f-ST(f)   e  k-ST(k)  =  C2  CT2     k"2/3  (25) 


where     S  (f)     and     ST(f)     are  spectral  density  values  with  units  of 

2         2 

(m/s)  /Hz  and  °C  /Hz  respectively.  The  above  equations  are  rewritten 

? 
in  the  following  forms  to  obtain  e  and  CT   from  values  of  f,  U, 


Su(f),  and  ST(f) 


wh 


e  =  (£ f  •  Su(f)}3/2  (26) 


CT2  =  {£ f  •  ST(f)}3/2  (27) 


ere  k  =  27rf/U  and  C,  =  .5  and  C«  =  .25  ,  empirically, 


C.  ANALYSES  OF  PAIRED  TEMPERATURE  SENSOR 

2 
Cj   values  were  also  estimated  by  using  paired  sensors  for  variance 

analyses.  This  analysis  was  based  on  the  expression  for  the  tempera- 
ture structure  function,  equation  (2). 

The  analysis  is  one  in  which  temperature  differences  between  paired 

resistance  wires  separated  by  a  distance  r  is  subjected  to  variance 

2 

analyses  as  indicated  by  the  term  [T1 (x)-T' (x+r)]  on  equation  (2). 
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Voltages  corresponding  to  the  temperature  differences,  AT,  as 
measured  by  the  paired  sensors  were  recorded  on  analog  magnetic  tape. 
During  experimental  periods,  calibration  voltages  were  also  recorded 
on  the  data  tape  for  reconstruction  of  the  signal.  This  calibration 
signal  was  played  back  at  the  time  of  data  analysis  to  provide  a 
reliable  check  of  the  recording  system  performance  in  both  record 
and  playback  modes. 

The  variance  or  mean  square  of  the  temperature  differences  required 

2 
to  compute  CT   was  obtained  by  analyzing  the  analog  voltages  with  a 

mean  square  voltmeter  (Thermo  Systems  Incorporated,  model  1060)  using 
a  30  second  time  constant  for  the  averaging.  The  output  of  the  mean 
square  voltmeter  as  well  as  a  record  of  the  original  signal  were  recorded 
on  standard  strip  charts.  A  strip  chart  section  appears  in  Figure  17. 
The  record  of  the  original  instantaneous  signal  provides  an  important 
means  of  spotting  extraneous  noises.  The  rms  voltmeter-strip  chart  re- 
corder system  was  also  calibrated  using  known  input  voltages.  At  least 
one  input  calibration  voltage  was  larger  than  the  largest  expected  signal 
on  the  data  record  and  likewise  one  input  calibration  voltage  was  lower 
than  the  smallest  expected  signal.  This  procedure  minimized  any  non- 
linearities  in  either  the  rms  voltmeter  of  the  strip  chart  recorder. 

Five  minute  averaged  values  of  the  rms  voltages  were  abstracted  from 
the  strip  chart  and  these  values  are  utilized  in  a  program  for  the  IBM 

360  computer.  This  program  converts  the  mean  square  voltages  values  to 

2 
temperature  difference  variances  values,  from  which  CT   was  computed. 

A  final  check  on  the  validity  of  the  measured  data  was  made  by  perform- 
ing spectral  analyses  on  representative  data  samples  from  each  experi- 
mental run.  The  spectrum  analyzer,  described  in  section  B,  was  used  to 
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Figure  17.  Sample  strip  chart  section, 
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obtain  the  spectra  which  were  subsequently  plotted  and  kept  as  a 
permanent  record.  Any  harmonic  noise  appeared  as  a  spike  in  an  other- 
wise reasonably  smooth  curve.  It  should  be  emphasized  that  continuous 
calibration  and  data  checking  procedure  were  essential  steps  through- 
out the  analyses. 

D.  HOT  WIRE  CALIBRATION 

In-situ  calibrations  of  the  velocity  sensors  were  accomplished 
every   10  minutes  during  the  experiment.  In  this  procedure  recordings 
were  made  of  both  the  cup  anemometer  wind  speed  and  the  corresponding 
hot  wire  voltage  output.  The  sensor  wind  speed  is  given  by 


2  =  a  U1/2  +  b  (28) 


V  =  a  U 

where     V     is  the  hot-wire  voltage  output  and     U     is  the  mean  wind  speed 
for  any  given  level.     The  constants,  a  and  b,  are  the  calibration  curve 
slope  and  intercept  respectively.     The  constants,  a  and  b,  are  derived 
from  the  in-situ  calibration  curve. 

To  convert  the  voltage  PSD  levels  to  velocity  units   requires  a 
calibration  factor  given  by 

< 

U'    =  Cv'  (29) 

where     C     is  the  calibration  factor  in  m/s  /volt,     v'     is  the  voltage 
fluctuation,  and     0"     the  velocity  fluctuation.     Differentiating  equa- 
tion (28)  yields 


"u"   =  (4vU1/2/a)   v'  (30) 
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which  when  substituted  into  equation  (29)  yields 


C  =  4vU1/2/a  (31) 


An  example  in-situ  calibration  curve  is  shown  in  Figure  18.  The 
slope  of  the  curve  yields  the  value  for  a  in  equation  (31). 
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Figure  18.     In-situ  calibration  plot, 
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V.   RESULTS 

A  comparison  of  the  observed  overwater  results  and  those  obtained 
by  Wyngaard,  et  al .  (1971)  appear  in  Figures  19,  20,  respectively. 
In  each  case  dots  represent  an  individual  data  point.  For  the  overwater 
results  mean  values  over  Ri  intervals  of  0.25  were  determined,  plotted 
in  the  center  of  the  interval,  and  appear  as  small  dots  within  circles. 
Error  bars  delineate  one  standard  deviation  from  the  mean  within  each 
interval.  The  number  at  the  top  of  each  error  bar  is  the  number  of 
data  points  used  in  defining  the  mean  in  the  interval. 

For  the  stable  and  unstable  regime  with  0.5  >  Ri  >-0.5  there  is 
some  correlation  with  the  results  of  Wyngaard,  et  al.  However,  for 
Ri  <-0.5  there  is  little  agreement  with  the  difference  in  results  on 
the  order  of  one  order  of  magnitude.  There  is  some  agreement  with  the 
trend  in  the  sense  that  it  does  not  change  significantly  with 
instability. 

Many  of  the  Richardson  numbers  were  not  within  acceptable  ranges, 
for  positive,  and  negative  Ri  values.  Results  are  not  included  for 
positive  Ri  values  above  1.0  and  negative  values  lower  than  -3.0.  A 
predicted  critical  positive  Ri  value  can  be  derived  from  equation  (19), 
which  gives  a  value  of  0.21. 

Scatter  in  the  observed  results  could  be  attributed  to  the  scatter 

2 
in  both  the  measured  CT   and  36/ 3z  values.  Deviation  of  tempera- 

2 
ture  spectra  from  the  -5/3  slope  caused  uncertainty  in  the  CT   esti- 
mates. The  30/82  values  are  dependent  on  profile  estimates  whose  de- 
termination was  at  times  subjective  as  discussed  in  conjunction  with 
Figure  21. 
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Figure  19.  Wyngaard,  et  al.  (1971)  prediction' curve. 


50 


Ct   /z  %  idd/dZ)2 


b     o 


O 

en       "_l       io 


01      _,       ro 


CO- 


01      O 

_l L_ 


1 
N3 

cn" 


»-s — G>- 


-©. 


CO 


10- 


-<2> H 


1 
01 


-© 


••  • 


■O-^ — 10 


•  ••       • 

-© H 


T-Z-& 


cn    • 


.  1 

01 


-»0-i ■— 1  tl 


•     • 


•  •  • 


3 )  cn 


t<& — *• 


*- : ^y 


-e- 


HS 


•  •  • 


>.  ». 


cn- 


•       • 


-^8> ; l£ 

I 0; ico 


_i«  • 


Figure  20.     Observed  prediction  curve, 
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Friction  velocity  values,  U*  ,  were  calculated  in  two  ways.  In 
the  first  method  (Method  A)  profiles  of  log  (z)  versus  mean  wind,  U  , 
were  plotted  to  obtain  gradient  estimate.  An  example  profile  is 
shown. in  Figure  10.  After  the  profiles  were  drawn,  5  and  10  meter 
level  wind  speeds  were  estimated  from  the  profiles  for  each  10  minute 
period  and  these  values  were  used  to  solve  equation  (14)  for  U*  . 
These  values  appear  in  column  A  of  Table  I. 

Friction  velocity  values  were  also  calculated  from  fluctuation 
data  (Method  B).  First,  equation  (24)  was  solved  to  obtain  a  value 
for  the  turbulent  kinetic  energy  dissipation  rate,  e  ,  at  each  level. 
Then  the  neutral  prediction,  equation  (16),  was  used  to  obtain  a  U* 
at  each  level.  Ideally  U*  at  all  levels  would  be  equal  for  any 
given  period  since  U*  is  assumed  to  be  constant  in  the  surface  layer, 
The  values,  with  about  10%  exceptions,  were  within  acceptable  limits 
of  each  other.  They  were  then  averaged  to  obtain  representative 
values  for  U*  .  Results  from  these  computations  appear  in  Table  I 
column  B. 

U*  values  were  also  calculated  using  wind  speed  dependent  C,n 
values  suggested  by  Cardone  (1969)  and  a  value  of  U  taken  from  pro- 
files, at  a  height  of  10  meters,  using 


U*  =  "10^  (32) 


Cardone' s  C,Q  representation  was  for  the  marine  boundary  layer. 
These  values  are  shown  in  Table  I  column  C. 

The  results  are  similar  to  those  reached  by  Atkinson  (1976),  i.e., 
showing  some  disagreement  between  U*  values  obtained  by  the  two 
methods,  with  the  fluctuation  data  method  more  closely  agreeing  with 
Cardone' s  results. 
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TABLE   I 
Comparison  of     U*     values. 


u. 


A 

B 

c 

.29 

.20 

.18 

.29 

.19 

.18 

.27 

.19 

.17 

.21 

.21 

.15 

.30 

.21 

.14 

.26 

'.21 

.14 

.25 

.22 

.14 

.21 

.26 

.17 

.19 

.24 

.15 

.19 

.24 

.15 

.08 

18 

.15 

.15 

.16 

.15 

.26 

.20 

.15 

.18 

.17 

.12 

A-      k(U2-Ui)/ln(z2/2l) 


METHOD:      B-      (*kz) 


V, 


Ui0(c10)/2 
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_3 
The  normalized  dissipation  rate  KcZU*    in  relation  to  Z/L  was 

tested  in  the  same  manner  as  Garratt  (1972).  The  results  of  Garratt 
(1972)  and  those  obtained  in  this  study  are  shown  in  Figures  22  and  23, 
respectively.  The  curve  represents  the  function  <j>  (Z/L) -Z/L  .  Similar 
scatter  was  also  observed  in  an  investigation  by  Busch  and  Panofsky 
(1968). 

An  examination  was  also  made  of  the  variation  of  the  observed  dissi- 
pation rate  of  turbulent  kinetic  energy,  e  ,  with  height  for  the  differ- 
ent stability  conditions.  Results  of  this  appear  in  Figures  24  and  25 
for  stable  and  unstable  conditions,  respectively.  The  results  from 
both  stability  regimes  are  compared  to  the  -1  properties  for  neutral 
stability.  The  results  show  that  for  all  stability  conditions  a  slope 
corresponding  to  near  neutral  conditions  was  observed. 
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30  APRIL  1976 


Figure  21.     Erroneous  profile. 
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Figure  22.     Garratt's  results 
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Figure  23.     Observed  results. 
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Figure  24.     e     versus  log  Z, 
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Figure  25.     e     versus  log  Z. 
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VI.  CONCLUSIONS 

2 
Correlation  of  spectrally  derived  CT   results  with  the  stability 

parameter,  Ri ,  was  poor  for  the  data  obtained  in  this  study  of  the 
marine  environment.  Verification  evaluations  for  friction  velocity 
results  from  mean  (profile)  and  turbulent  (dissipation)  data  were  satis- 
factory. Correlation  of  normalized  e  with  the  stability  parameter 
(Z/L)  showed  the  same  scatter  as  results  of  others.  The  variation  of 
e  with  height  exhibited  results  which  support  the  near  neutral  pre- 
dictions for  most  cases. 

It  is  concluded  that  the  turbulent  field  in  the  marine  boundary 
layer  is  subjected  to  several  anomalous  effects  which  can  cause  the 
disruption  of  the  inertial  subrange.  Thus,  empirically,  expressions 
for  describing  overland  transmission  characteristics  will  undoubtedly 
have  to  be  altered  to  provide  operational  determinations  for  utiliza- 
tion in  the  marine  environment. 
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1974 


TABLE   II. 

2 
Cy     and  Ri   Results 

• 

Period 

Level 
(M) 

de/dz 
°C/cm 

c  2 

T     -3 
x'lO  J 

Ri 

3-27-1841 

4.19 

6.60 

7.64 

13.89 

0.049 
0.031 
0.027 
0.015 

1.59 
0.94 
1.57 

0.362 
0.500 
0.500 
1.000 

3-27-1902 

4.19 

6.60 

7.64 

13.89 

0.047 
0.030 
0.026 
0.014 

1.97 
0.61 
2.02 

0.249 
0.392 
0.454 
0.750 

3-27-1928 

4.19 

6.60 

7.64 

13.89 

0.047 
0.030 
0.026 
0.014 

5.40 
1.37 
2.89 
1.20 

0.230 
0.362 
0.419 
0.750 

3-27-2117 

4.19 

6.60 

7.64 

13.89 

0.019 
0.012 
0.012 
0.006 

3.45 
6.89 

1.42 

0.073 
0.115 
0.113 
0.242 

3-28-1902 

4.19 

6.60 

7.64 

13.89 

0.046 
0.029 
0.025 
0.014 

16.56 

12.07 
1.46 

1.000 
1.500 
1.500 
1.500 

3-28-2059 

4.19 

6.60 

7.64 

13.89 

0.054 
0.034 
0.030 
0.016 

56.73 

14.21 
19.25 

0.388 
0.500 
0.500 
1.000 

8-13-1700 

4.19 

6.60 

7.64 

13.89 

-0.036 
-0.023 
-0.020 
-0.011 

4.33 
0.42 
0.53 
0.23 

-0.538 
-0.847 
-0.981 
-1.783 

8-13-1721 

4.19 

6.60 

7.64 

13.89 

-0.036 
-0.023 
-0.020 
-0.011 

3.76 
0.42 
0.67 
0.27 

-0.665 
-1.047 
-1.213 
-2.204 
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1974 


Period 

Level 
(M) 

de/dz 
°C/cm 

c2 

xT10"3 

Ri 

8-13-1742 

4.19 

6.60 

7.64 

13.89 

-0.025 
-0.016 
-0.013 
-0.007 

4.25 
0.45 
0.63 
0.21 

-0.451 
-0.711 
-0.823 
-1.496 

9-18-1643 

4.19 

6.60 

7.64 

13.89 

0.031 
0.020 
0.017 
0.009 

0.29 
0.08 

0.11 

0.308 
0.486 
0.500 
1.000 

9-19-1735 

4.19 

6.60 

7.64 

13.89 

-0.005 
-0.003 
-0.003 
-0.001 

0.47 
1.18 

0.56 

-0.082 
-0.128 
-0.149 
-0.270 

9-19-1846 

4.19 

6.60 

7.64 

13.89 

-0.005 
-0.003 
-0.003 
-0.001 

5.07 
10.27 

1.10 

-0.100 
-0.157 
-0.181 
-0.330 

9-19-1900 

4.19 

6.60 

7.64 

13.89 

-0.010 
-0.006 
-0.005 
-0.003 

0.59 
0.87 

0.91 

-0.607 
-0.957 
-1.107 
-2.013 

9-19-1929 

4.19 

6.60 

7.64 

13.89 

-0.005 
-0.003 
-0.003 
-0.001 

0.73 
1.17 

1.14 

-0.099 
-0.156 
-0.180 
-0.328 

9-19-2011 

4.19 

6.60 

7.64 

13.89 

-0.010 
-0.006 
-0.005 
-0.003 

0.50 
0.81 

0.55 

0.051 
0.080 
0.093 
0.169 

9-19-2115 

4.19 

6.60 

7.64 

13.89 

-0.172 
-0.109 
-0.094 
-0.052 

0.35 
0.69 

0.77 

-3.000 
-3.000 
-3.000 
-3.000 
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1974 


Period 

Level 
(M) 

de/dz 
°C/cm 

c2 

T     -3 
x'lO  J 

Ri 

9-19-2158 

4.19 

6.60 

7.64 

13.89 

-0.010 
-0.006 
-0.005 
-0.003 

0.46 
1.12 

1.10 

-0.202 
-0.319 
-0.369 
-0.671 

9-19-2243 

4.19 

6.60 

7.64 

13.89 

-0.178 
-0.113 
-0.098 
-0.054 

0.32 
0.78 

1.08 

-3.000 
-3.000 
-3.000 
-3.000 

9-20-0014 

4.19 

6.60 

7.64 

13.89 

-0.011 
-0.007 
-0.006 
-0.003 

0.66 
2.12 

1.0 

-0.461 
-0.727 
-0.841 
-1.530 

9-20-0056 

4.19 

6.60 

7.64 

13.89 

-0.013 
-0.008 
-0.007 
-0.004 

0.56 
1.88 

-0.408 
-0.643 
-0.744 
-1.353 

9-20-0140 

4.19 

6.60 

7.64 

13.89 

-0.013 
-0.008 
-0.007 
-0.004 

0.77 
1.43 

0.89 

-0.546 
-0.860 
-0.995 
-1.809 

9-20-0530 

4.19 

6.60 

7.64 

13.89 

-0.003 
-0.002 
-0.002 
-0.001 

0.76 
1.77 

1.68 

-0.085 
-0.134 
-0.155 
-0.282 

9-20-0552 

4.19 

6.60 

7.64 

13.89 

-0.006 
-0.004 
-0.004 
-0.002 

0.91 
0.99 

-0.952 
-1.500 
-1.736 
-3.000 

11-21-1800 

4.19 

6.60 

7.64 

13.89 

-1.626 
-1.032 
-0.892 
-0.490 

3.74 
8.85 
2.02 
5.93 

-3.000 
-3.000 
-3.000 
-3.000 
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1974 


1975 


Period 

Level 
(M) 

d0/dz 
°C/cm 

c  2 
x  10 

Ri 

11-21-1821 

4.19 

6.60 

7.64 

13.89 

0.008 
0.005 
0.004 
0.002 

3.39 
1.97 
1.97 
2.24 

0.027 
0.043 
0.050 
0.091 

11-21-1842 

4.19 

6.60 

7.64 

13.89 

0.006 
0.004 
0.004 
0.002 

3.49 
2.17 
2.11 
9.17 

0.033 
0.051 
0.060 
0.108 

11-21-1903 

4.19 

6.60 

7.64 

13.89 

0.003 
0.002 
0.002 
0.001 

5.65 
3.26 
3.00 
1.24 

0.009 
0.014 
0.017 
0.030 

11-21-2027 

4.19 

6.60 

7.64 

13.89 

0.002 
0.001 
0.001 
0.000 

7.75 
1.42 
3.61 
0.89 

-0.001 
-0.001 
-0.001 
-0.002 

11-21-2048 

4.19 

6.60 

7.64 

13.89 

0.002 
0.001 
0.001 
0.000 

8.39 
9.75 
3.55 
1.36 

0.006 

0.010 

-   0.012 

0.021 

11-21-2109 

4.19 

6.60 

7.64 

13.89 

-0.000 
-0.000 
-0.000 
-0.000 

9.71 
7.22 
4.11 
1.19 

-0.025 
-0.039 
-0.046 
-0.083 

11-21-2130 

4.19 

6.60 

7.64 

13.89 

0.0 
0.0 
0.0 
0.0 

11.04 
7.97 
5.90 
1.72 

-0.037 
-0.059 
-0.068 
-0.123 

3-27-1503 

4.19 

6.60 

7.64 

13.89 

0.006 
0.004 
0.004 
0.002 

10.57 
6.14 
3.43 

0.004 
0.007 
0.008 
0.015 
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1975 


Period 

Level 
(M) 

de/dz 
°C/cm 

C  2 
x  10 

Ri 

3-27-1519 

4.19 

6.60 

7.64 

13.89 

0.006 
0.004 
0.004 
0.002 

9.66 
5.90 
1.65 

,     0.004 
0.006 
0.007 
0.012 

3-27-1539 

4.19 

6.60 

7.64 

13.89 

-0.002 
-0.001 
-0.001 
-0.000 

10.04 
5.86 
2.58 

-0.059 
-0.092 
-0.107 
-0.195 

3-27-1600 

4.19 

6.60 

7.64 

13.89 

-0.005 
-0.003 
-0.003 
-0.001 

10.71 
3.42 
2.65 

-0.032 
-0.051 
-0.059 
-0.107 

3-27-1644 

4.19 

6.60 

7.64 

13.89 

-0.002 
-0.001 
-0.001 
-0.000 

11.25 
7.46 
3.01 

-0.014 
-0.023 
-0.026 
-0.048 

3-27-1656 

4.19 

6.60 

7.64 

13.89 

0.0 
0.0 
0.0 
0.0 

11.31 
5.74 
2.08 

-0.013 
-0.020 
-0.024 
-0.043 

3-27-2020 

4.19 

6.60 

7.64 

13.89 

0.041 
0.026 
0.022 
0.012 

2.76 

0.016 
0.025 
0.029 
0.053 

3-27-2045 

4.19 

6.60 

7.64 

13.89 

0.044 
0.028 
0.024 
0.013 

5.17 
0.20 

0.025 
0.039 
0.045 
0.082 

3-27-2106 

4.19 

6.60 

7.64 

13.89 

0.015 
0.009 
0.008 
0.004 

5.45 
2.01 

-0.046 
-0.072 
-0.083 
-0.151 
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TABLE   III 

April   Cru' 

2 
ise  Cj     and  Ri 

Date  Level    1 

Results 

Time 

80/8Z 

2       3 
CT  xlO^ 

Ri 

(°C/m) 

4-27-1519 

.135 

4.439 

-0.68 

1529 

.150 

3.288 

-1.15 

1544 

.160 

2.877 

-0.88 

1554 

.160 

2.466 

-0.73 

1625 

.160 

2.634 

-1.42 

1635 

.175 

2.404 

-2.04 

4-28-1433 

.050 

5.178 

-0.75 

1443 

.075 

4.439 

-0.93 

1458 

.050 

7.398 

-0.60 

1517 

.085 

4.809 

-0.80 

1531 

.085 

5.918 

-1.88 

1540 

.050 

5.918 

-20.01 

1613 

.075 

5.548 

-0.94 

1618 

.075 

5.548 

-0.90 

1623 

.105 

4.439 

-1.84 

1633 

.065 

5.918 

-1.98 

1703 

.035 

2.877 

-0.37 

1720 

.125 

3.288 

5.95 

1730 

.925 

3.699 

0.70 

1740 

.020 

3.699 

-0.47 

1909 

.040 

2.774 

1.07 

1919 

.040 

2.959 

4.45 

1929 

.030 

3.144 

5.01 

4-29-1415 

.025 

1.099 

3.34 

1425 

.020 

3.625 

0.50 

2035 

.040 

3.625 

-0.99 

2050 

,040 

3.625 

-1.48 

4-30-1315 

.040 

4.285 

-0.71 

1325 

.040 

5.933 

-1.80 

1335 

.050 

5.273 

-0.61 
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Data  Level   2 

Time 

dQ/dz 

2       3 
CT  xlO"3 

Ri 

(°C/m) 

4-27-1519 

.095 

4.279 

-1.07 

1529 

.105 

2.445 

-1.81 

1544 

.10 

2.649 

-1.38 

1554 

.105 

2.852 

-1.15 

1625 

.105 

1.630 

-2.24 

1635 

.11 

1.936 

-3.22 

4-28-1433 

_ 

_ 

-1.18 

1443 

.055 

3.260 

-1.46 

1458 

- 

- 

-0.95 

1517 

.060 

4.482 

-1.26 

1531 

.060 

3.871 

-2.96 

1540 

.035 

3.871 

-31.52 

1613 

.050 

3.667 

-1.48 

1618 

.055 

4.279 

-1.42 

1623 

.065 

3.464 

-2.89 

1633 

.015 

4.279 

-3.12 

1703 

.030 

2.649 

-0.58 

1720 

.955 

3.056 

9.37 

1730 

.001 

3.260 

1.11 

1740 

.025 

3.056 

-0.74 

1909 

.01 

2.751 

1.68 

1919 

.015 

2.037 

7.01 

1929 

.005 

- 

7.89 

4-29-1415 

.03 

4.265 

5.26 

1425 

.005 

4.436 

0.79 

2035 

.035 

1.706 

-1.57 

2050 

- 

- 

-2.34 

4-30-1315 

.020 

2.730 

-1.12 

1325 

.020 

2.218 

-2.83 

1335 

.050 

3.071 

-0.96 
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Data 

Level   3 

Time 

86/3z 

2       3 
CT  xlOJ 

Ri 

(°C/m) 

4-28-1340 

.045 

3.129 

-1.69 

1350 

.045 

2.869 

-1.37 

1400 

.040 

2.608 

-1.24 

1410 

.040 

2.869 

-0.91 

1423 

.035 

2.347 

-1.21 

1803 

.035 

0.652 

-0.68 

1814 

.035 

0.652 

-0.68 

1827 

.015 

0.782 

-0.61 

1836 

.01 

0.782 

0.27 

1847 

.01 

3.912 

0.34 

1858 

.01 

9.779 

0.34 

4-29-1700 

.065 

2.539 

-7.12 

1711 

.08 

2.537 

-4.7 

1722 

.075 

3.173 

-12.6 

1732 

.08 

3.173 

-5.98 

1744 

.065 

2.539 

-5.98 

4-30-1125 

.05 

3.491 

-2.92 

1137 

.015 

2.856 

0.43 

1232 

.01 

1.904 

0.0 

1242 

.015 

1.523 

0.0 

1252 

.01 

1.142 

0.65 

1302 

.015 

1.269 

0.0 
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Data  Level  4 

Time 

86/8Z 

CT2xl03 

Ri 

(°C/m) 

4-27-1519 

.075 

2.658 

-2.24 

1529 

.075 

1.772 

-3.82 

1544 

.080 

1.772 

-2.90 

1554 

.085 

1.772 

-2.41 

1625 

.080 

1.329 

-4.71 

1635 

.085 

1.476 

-6.77 

4-28-1433 

.015 

4.958 

-2.49 

1443 

.035 

2.953 

-3.07 

1458 

.020 

6.325 

-1.99 

1517 

.040 

4.616 

-2.65 

1531 

.040 

2.067 

-6.22 

1540 

.025 

2.658 

-66.33 

1613 

.030 

2.658 

-3.11 

1618 

.035 

2.953 

-2.99 

1623 

.055 

2.067 

-6.09 

1633 

.045 

- 

-6.57 

1703 

.015 

1.772 

-1.23 

1720 

.075 

2.067 

19.71 

1730 

.020 

2.067 

2.33 

1740 

.010 

1.772 

-1.55 

1909 

.040 

2.137 

3.54 

1919 

.025 

1.476 

14.76 

1929 

.025 

1.476 

16.61 

4-29-1415 

.005 

2.531 

11.06 

1425 

.015 

2.847 

1.65 

2035 

.005 

- 

-3.29 

2050 

.005 

- 

-4.92 

4-30-1315 

.015 

_ 

-2.35 

1325 

.015 

- 

-5.96 

1335 

.005 

- 

-2.03 
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TABLE  IV 

3 
April  Cruise  e  Results  (e  x  10  ) 


Time       Level  1       Level  2       Level  4 


4-28-1330 

3.61 

4.92 

1.62 

1339 

3.97 

2.90 

1.66 

1358 

3.58 

3.67 

1.30 

1757 

4.37 

4.55 

1.93 

1815 

5.83 

4.66 

1.56 

1833 

4.68 

4.95 

1.88 

1851 

4.33 

3.94 

1.87 

4-29-1652 

12.4 

4.00 

6.68 

1710 

8.41 

2.33 

5.83 

1728 

10.4 

3.04 

4.17 

1916 

6.92 

11.6 

2.44 

1934 

5.04 

1.30 

1.80 

2048 

- 

2.21 

4.35 

2203 

17.8 

3.12 

1.15 

2221 

19.4 

3.43 

1.78 

4-30-1013 

4.87 

3.43 

0.84 

1031 

4.82 

0.97 

0.88 

1232 

7.40 

2.09 

1.93 

1250 

3.67 

1.37 

1.37 
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